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Abstract. We present here a computer model of the dis- 
tribution of the luminosity, produced by old isolated neu- 
tron stars and black holes accreting from the interstellar 
medium in the Galaxy in R—Z plane. We show that the lu- 
minosity distributions have a toroidal structure, with the 
maximum at w 5 — 6kpc. 
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1. Introduction 

Old isolated neutron stars (NS) and black holes (BH) form 
a big populations of the galactic sources (about 10 8 -10 9 
objects in the Galaxy), but most of them are unobserved 
today. Less than 10 3 of young NS appear as radiopulsars, 
and no isolated BH are observed at the present moment 
(probably, some of them were detected, for example, in 
the ROSAT survey, but no one is identified). As in this 
article we will speak only about isolated compact objects, 
we will not use that adjective: "isolated" for NS and BH. 

During the last several years, the spatial distribution 
and other properties of NS became of great interest, be- 
cause as it has been proposed, NS can be observed by 
the ROSAT sattelite in soft X-rays due to accretion from 
the interstellar medium (ISM) (see, for example, Treves & 
Colpi 1991), and several sources of this type have been ob- 
served (Walter et al. 1996). BH also can appear as such X- 
ray sources (Heckler & Kolb 1996) with some differences 
in spectrum and temporal behaivour (absence of pulsa- 
tions, for example). That's why here we try to obtain a 
picture of the distribution of the accretion luminosity of 
these sources. 

Fast rotation and/or strong magnetic field can prevent 
accretion onto the surface of the NS. In this case the X- 
ray luminosity will be very low (except the case of the 
transient source formation due to the envelope around the 
NS) (see Popov 1994 and Lipunov & Popov 1995). Here we 
consider only accreting NS. Most of NS are on the stage 



of accretion, because their magneto-rotational evolution 
usually finishes at this stage approximately 10 8 years after 
their birth. The NS properties (periods etc) on the stage 
of accretion depend upon the magnetic field decay (see 
Koncnkov & Popov 1997). BH, of course, can appear only 
as Accretors. 

In the articles of Gurevich et al. (1993), Postnov & 
Prokhorov (1993, 1994) it was shown that NS in the 
Galaxy form a toroidal structure. The distribution of 
the ISM (see, for example, Bochkarev 1993) also has 
the toroidal structure. The maxima of both distributions 
roughly coincide. 

Therefore, most part of NS (and, as one can say, BH) is 
located in the dence regions of the ISM. Thus the accretion 
luminosity in these regions should be higher. The results 
of computer simulations of this situation are presented in 
this paper. 

The trajectories of NS and BH were computed directly 
for specified initial velocity distribution, the Galaxy gravi- 
tational potential and the distribution of the ISM density. 
Preliminary results of such computations for NS for 8- 
function and maxwellian velocity distributions were pre- 
sented in Popov & Prokhorov (1996, paper I). 

In the section 2 we briefly describe our model. In the 
third section the results and a short discussion are pre- 
sented. And finaly we conclucd our results in the last sec- 
tion. 

2. The Model 

We solved numerically the system of differential equations 
of motions in the given potential. 

We used the Galactic potential, taken in the form 
(Paczynski 1990): 



^(R, Z) = GM t / (r 2 + [a, + {Z 2 + 6 2 ) 1/2 ] 2 ) 
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with a quasi-spherical halo with the density distribution 
in the form: 
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Disk 

Bulge 

Halo 



a D =0 
<2b=3.7 kpc 



M D = 1.12 • 1O 1O M 
Mb = 8.07 • 10 10 M Q 
M = 5.0 • 1O 1O M 



Here R and Z are cylindrical coordinates, d- radial 
distance in the quasi-spherical halo. The parameters of the 
potential are given in the table, po is determined through 
the halo mass. M . 

b D =277 pc 
fe B =200 pc 
d =277 pc 

The density in our model was constant in time. Lo- 
cal density was calculated using data and formulae from 
Bochkarev (1993) and Zane et al. (1995). Here R and Z 
are cylindrical coordinates (it is assumed, that the Galaxy 
has cylindrical symmetry), n - total concentration, nni 
and nn 2 are concentrations of the neutral and molecu- 
lar hydrogene, no(R), U2(R) and n 3 (i?) are the values of 
concentrations for Z = 0. 



n(R, Z) = n HI + 2 • n H , 2 
-Z 2 



Dan ! 



nn 2 = n 2 (R) ■ exp 
If Okpc <R< 3Akpc, then 

nni = no(R) ■ exp 



2 • (70pc) 2 



-Z 2 



2 • (UOpc ■ R/2kpcf 

For Okpc < R < 2kpc no(R) was assumed to be con- 
stant for all R: 

no(R < 2kpc) = n(R = 2kpc) 

Of course, it is not accurate for small R, so for the 
very central part of the Galaxy our results arc only a 
rough estimation (see Zane et al. (1996) for detailed cal- 
culation of NS emiision from the Galactic center region). 
If 3Akpc <R< 8.5fcpc, then 



nni — 0.345 • exp 



0.107- exp 



-Z 2 



2 ■ (212pc) 2 
-Z 2 



+ 



0.064 • exp 
If 8.5 < R < Wkpc, then 

nni = n 3 (R) ■ exp 



2 ■ (530pc) 2 

Z " 



+ 



403pc 



-Z 2 



2 ■ (530pc • R/S.bkpc) 2 

no(R), ri2 and ^(R) were taken from Bochkarev (1993). 

The density distribution in the R-Z plane used in our 
computations is shown in the figure 1. 

In our model we assumed, that the birthrate of NS 
and BH is proportional to the square of the local density. 
Stars were born in the Galactic plane (Z=0) with circular 
velocities plus additional isotropic kick velocities. 
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Fig. 1. The density distribution in units particle per cubic 
centimeter in the R-Z plane. 



For the kick velocity we used the formula from Lipunov 
et al. (1996). It was constructed as an analitical approxi- 
mation of the three-dimensial velocity distribution of ra- 
diopulsars from Lyne & Lorimer (1994). 



„0.19 



fLL(V) CX 



( 1 + 2-6.72)1/2' 



here V- spatial velocity of the compact object, V c har - 
characteristic velocity, x = V/V c har, III- the probability 
(see the detail describtion of the analitical approximation 
in Lipunov et al. (1996)). This formula reproduces the ob- 
served distribution with the mean velocity 350 km/s for 
V r c /iar=400 km/s. This velocity distribution seems more 
luckely, than the 5-function and Maxwellian distributions, 
which we used in Paper I. Kick velocities were taken equal 
for NS and BH (this is a "zero-hypothesis", as we don't 
have any exact information on this matter). But it is pos- 
sible, that BH have lower kick velocities because of their 
high masses. One of the reasons to make computations for 
V c har— 200 km/s was to explore this situation. 

For each star we computeded the exact trajectory and 
the accretion luminosity. The accretion luminosity was cal- 
culated using Bondi formula: 



L 



( gmm\ 
\ Rub J 



■ f (GM) 2 p(R,Z) 



Sound velocity, V s , was taken to be 10 km/s every- 
where. We used equal masses M^s — 1-4M Q for NS 
and Mbh = 10M Q for BH. Density is determined as: 
p = n ■ niH, where m#- the mass of the hydrogen atom. 
Radii, Rub, where the energy is libirated, were assumed 
to be equal to 10 km for NS and 90 km (i.e. 3 • R g , 
R g = 2GM/c 2 ) for BH. Statistics was collected on the 
grid with the cell size 100 pc in R-direction and 10 pc in 
Z-direction (centered at R=50 pc, Z=5 pc and so on). The 
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luminosity is shown on the graphs in ergs per second per 
cubic parsec. 

For the normalization of our results we used an as- 
sumption, that there are Nns = 10 9 NS and Nbh = 10 8 
BH in the calculated volume of the Galaxy. For Salpeter 
mass function with a=2.35 the ratio of NS to BH is about 
10 if NS are formed from the stars with masses between 
10M Q and w 45 - 50M Q , and BH from the stars with 
masses higher than « 45 — 50M Q . Motch et al. (1997) ar- 
gued, that N ns — 10 9 can be ruled out, and N NS = 10 8 
is a more probable value, but for the caclculations of the 
distribution the total number is not so improtant, and 
for other numbers of compact objects the results (i.e. the 
value of the luminosity) can be easily scaled. 
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Fig. 3. The accretion luminosity distribution in the R-Z plane 
for neutron stars for characteristic kick velocity 400 km/s. 
The luminosity is shown in ergs per second per cubic parsec. 
N NS = 10 9 



3. Results and discussion 

On the figures 2-5 we represent the results for two charac- 
teristic values of the velocity distribution for NS and BH. 
On the graphs 2-5 the region 16000 pc x 800 pc for aver- 
aged accretion luminosity for the azimutale angle 0-180 
degrees is shown. We show only this region, because we 
are mainly interested in the part, where the maximum is 
situated, but while the computations NS and BH could 
move inside infinite region. The scale for R and Z axes 
are different in order to show clearly the structure in Z 
direction. Differences between the luminosity distribution 
for Z > and Z < demonstrate accuracy of the compu- 
tations (curves were not smoothed). 
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Fig. 2. The accretion luminosity distribution in R-Z plane 
for neutron stars for characteristic kick velocity 200 km/s. 
The luminosity is shown in ergs per second per cubic parsec. 
Nns = 10 9 

On the figure 6 the slice at Z=+5 pc (the first grid 
cells in the positive Z-direction) for the characteristic kick 
velocity V c har~ 200 km/s is shown: luminosity vs. radius. 
The figure is not symmetric. The right part corresponds 
to the azimutale angles 0-180 degrees, the left - 180-360 
degrees. Differences between the left and the right parts 
of the curve demonstrate how precise our method is. 
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R , pL- 
Fig. 4. The accretion luminosity distribution in the R-Z 
plane for black holes for characteristic kick velocity 200 km/s. 
The luminosity is shown in ergs per second per cubic parsec. 
Nbh = 10 s 

4. Concluding remarks 

As it is clearly seen from the figures, the distribution of 
the accretion luminosity in R-Z plane forms a toroidal 
structure with the maximum at approximatelly 5 kpc. 

As expected, for BH we have relatively higher lumi- 
nosity, than for NS, as BH have greater masses. But if the 
total number of BH is significantly lower than the num- 
ber of NS, their contribution to the luminosity can be less 
than the contribution of NS. The total accretion lumi- 
unosity of the Galaxy for N^s = 10 9 and Nbh = 10 8 is 
about 10 39 — 10 40 erg/s. For the characteristic velocity 200 
km/s for NS the maximum of the distribution is situated 
approximately at 5.0 kpc and for BH at 4.8 kpc. For NS 
with the characteristic velocity 400 km/s maximum is sit- 
uated at 5.4 kpc, and for BH at 5.0 kpc. This result is 
also expected because of the same reason (higher masses 
of the BH with the same value of kick velocity). 

The toroidal structure of the luminosity distribution 
of NS and BH is an interesting and important feature of 
the Galactic potential. As one can suppose, for low char- 
acteristic kick velocities and for BH we obtained higher 
luminosity. 
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Nbh = 10 s 
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Fig. 6. Slice at Z=+5 pc for NS for the characteristic kick 
velocity 200 km/s. Nns = 10 9 . The accretion luminosity is 
shown in ergs per second per cubic parsec. With the solid line 
we show the smoothed curve. 



As me made very general assumptions, we argue, that 
such a distribution is not unique only for our Galaxy, and 
all spiral galaxies can have such a distribution of the ac- 
cretion luminosity, associated with accreting NS and BH. 
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